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Short title: PHOTOSPHERIC MAGNETIC PRELUDE TO FLARES 
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Abstract. The observations of magnetic field variations as a signature of 
flaring activity is one of the main goal in solar physics. Some efforts in the past give 
apparently no unambiguous observations of changes. We observed that the scaling 
laws of the current helicity inside a given flaring active region change clearly and 
abruptly in correspondence with the eruption of big flares at the top of that active 
region. Comparison with numerical simulations of MHD equations, indicates that the 
change of scaling behavior in the current helicity, seems to be associated to a topological 
reorganization of the footpoint of the magnetic field loop, namely to dissipation of small 
scales structures in turbulence. It is evident that the possibility of forecasting in real 
time high energy flares, even if partially, has a wide practical interest to prevent the 
effects of big flares on Earth and its environment. 
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Solar flares are sudden, transient energy release above active regions of the Sun 
(Priest, 1982). As a consequence of random motion of the footpoints of the magnetic 
field in the photospheric convection (Parker, 1988), flares represent the dissipation 
at the many tangential discontinuities arising spontaneously in the magnetic fields of 
active regions. The magnetic energy is released in various form as thermal, soft and 
hard X-ray, accelerated particles etc. The observations of magnetic field variations, as a 
signature of flares above active regions, has been one of the main goals in solar physics, 
and some attempts for this has been made in the past (e.g. Hagyard et ai, 1999 and 
references therein). All efforts give apparently no unambiguous observations of changes. 
This is due to the fact that often investigations try to look for changes of the vector 
magnetic field as a whole. Recently, unambiguous observations of changing have been 
reported by Yurchyshyn et ai, 2000. The authors observed some typical changes of 
the scaling behavior of the current helicity calculated inside an active region of the 
photosphere, connected to the eruption of big flares above that active region. In the 
present paper we conjecture that the changes in the scaling behavior of the observed 
quantity is related to the occurrence of changes in the topology of the magnetic field at 
the footpoint of the loop. 

The occurrence of scaling of signed measures calculated from scalar fields f(x) 
which oscillate in sign can be studied through the following steps. First of all, introduce 
the signed measure 
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through a coarse-graining of non overlapping boxes Qi(r) of size r, covering the whole 
field defined on a region of size L. It has been observed (Ott et ai, 1992) that, for fields 
presenting self-similarity, this quantity displays well defined scaling laws. That is, in a 
range of scales r, the partition function x( r )> defined as 



where the sum is extended over all boxes occurring at a given scale r, follows a power- law 
behavior 



The scaling exponent k has been called cancellation exponent (Ott et ai, 1992) because 
it represents a quantitative measure of the scaling behavior of imbalance between 
negative and positive contributions in the measure. For a positive definite measure or a 
smooth field k — 0, while k = d/2 for a completely stochastic field in a d- dimensional 
space (for example a field of uncorrelated points with / = ±1, the sign being choosen 
randomly and independently for each point, with probability 1/2). As the cancellations 
between negative and positive part of the measure decreases toward smaller scales, we 
get k > 0, and this is the interesting situation. It is clear that the presence of structures 
in the field has an important effect on the cancellation exponent. For example, values of 
k < d/2, where d is the dimension of the space (in the present paper d — 2), indicate 
the presence of sign-persistent («. e. smooth) structures. 

Within turbulent flows, the value of the cancellation exponent can be related to the 
characteristic fractal dimension D of turbulent structures on all scales using a simple 




(2) 



Qi(r) 




(3) 
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geometrical argument (Sorriso-Valvo et ai, 2002). Let A be the typical correlation 
length of that structures, of the order of the Taylor microscale (see for example Frisch, 
1995), so that the field is smooth (correlated) in D dimensions with a cutoff scale A, 
and uncorrelated in the remaining d — D dimensions. If the field is homogeneous, the 
partition function (0) can be computed as (L/r) d times the integral over a generic 
box Q(r) of size r. The scaling of the latter can be estimated integrating over regular 
domains of size X d and considering separately the number of contributions coming 
from the correlated dimensions of the field and those from the uncorrelated ones. The 
integration of the field over the smooth dimensions will bring a contribution proportional 
to their area (r/X) D , while the uncorrelated dimensions will contribute as the integral of 
an uncorrelated field, that is proportional to the square root of their area (r / \p d ~ D >' 2 . 
Thus, when homogeneity is assumed, collecting all the contributions in (^) leads to 
scaling x( r ) ~ r -{d-D)/2 f or ^ e partition function, so that one can recover the simple 
relation 

K =(d-D)/2. (4) 

To get a quantitative measure of the change of the scaling of current helicity 
inside active regions, we used observations of the vector magnetic field obtained with 
the Solar Magnetic Field Telescope of the Beijing Astronomical Observatory (China). 
Measurements were recorded in the Fel 5324.19 A spectral line. The field of view is 
about 218" x 314", corresponding to 512 x 512 pixels on CCD. The magnetic field 
vector at the photosphere has been obtained through the measurements of the four 



7 



Stokes parameters, and the current density J z (x,y) has been calculated as a line 
integral of the transverse field vector over a closed contour of dimension 1.72" x 1.86" 
(cf. Yurchischin et al, 2000, for details). The current helicity H c = B ■ J (where B 
represents the magnetic field and J = V x B the current density) is a measure of 
small scales activity in magnetic turbulence. It indicates the degree of clockwise or 
anti-clockwise knotness of the current density. Let us consider a magnetogram of size L 
taken on the solar photosphere of an active region, and let T$±(x,y) the observed 
magnetic field perpendicular to the line of sight ((x,y) are the coordinate on the surface 
of the sun). Through this field we can measure the surrogate of current helicity, that is 
h c (x,y) = B z (x,y)J z (x,y) being J z (x,y) = [V x Bj_] ■ e z . Figure [I] shows the current 
helicity surrogate of the active region NOAA 7590 for five different times, during the 
flaring event. The presence of signed structures is evident, as well as their evolution 
with time. A signed measure can be defined from this quantity 



In figure [| we show, as example, the scaling behavior of x( r ) vs - r f° r a flaring active 
region (NOAA 7315) which started to flare on October 22, 1992. At larger scales we find 
x(r) ~ const., and this is due to the complete balance between positive and negative 
contributions. The same behavior does not appear at smaller scales, showing that the 
resolution of the images is not high enough to resolve the smallest structures. In the 
intermediate region of scales, the cancellation exponent is found to be k = 0.53 ± 0.09 
(Yurchischin et al, 2000). 




h c (x, y)dxdy . 



(5) 



Let us consider now what happens to the fractal dimension of current structures D 
as a function of time. To this aim, we take different consecutive magnetograms of the 
same active region, and for each magnetogram, we compute the value of k and then 
of D through relation (£|). Note that, since cancellations in the vertical photospheric 
magnetic field B z (x,y) itself have been found to be very small (Lawrence et al, 1993; 
Abramenko et al, 1998), with a cancellation exponent of the order of 10 -2 , cancellations 
of the current helicity are entirely due to the current structures. 



In figure [|, we report the time evolution of D superimposed to the flares occurred figure |3| 
in two active regions, namely NOAAs 7315, and 7590 (which flared on October 1, 1993). 
Quite surprisingly we observe that the fractal dimension D, starting from a given value 
(D < 1), becomes abruptly larger in correspondence with a sequence of big flares 
occurring at the top of the active region into the corona. The same behavior has 
been found for all calculations in all active regions we examined. The increase of the 
dimension of the structures may be the signature that dissipation has occured. In fact, 
annihilation is responsible for the smoothing of the small scales structures. 

As shown in figure ^|, a saturation of x( r ) is observed at large scales. At the 
smallest scales, the density of the measure must becomes smooth (no changes in sign 
are present) and we might thus found a saturation of x( r )- The fact that we do not 
find this saturation is an indirect evidence that elementary flux tubes are smaller than 
the instrumental resolution. The change towards larger D is then probabily due to the 
occurrence of dissipation of smaller and smaller flux tubes, that is magnetic energy is 
suddenly transferred towards small scales. This is the occurrence of an energy cascade 



towards smaller scales. 

Using high resolution numerical simulation of two-dimensional (d = 2) turbulent 
magnetohydrodynamic flows (Politano et ai, 1998; Sorriso-Valvo et ai, 2000; Sorriso- 
Valvo et ai, 2001), we can build up the signed measure for different fields. For example, 
since the geometry of the magnetic field B(x,y) = (B x ,B y ,0) is two-dimensional, the 
current J(x, y) = V x J = (0, 0, J z ) has only the z component, perpendicular to the 2-d 



simulation box, %. e. the plane (x,y). In Figure | we display the current field J(x,y) for Figure |] 

the numerical data, using ten snapshots in the statistically steady state, from t = 168 

up to t = 336 in non-linear times units, tnl- As can be seen, the presence of positive 

and negative structures is evident as in the case of the solar data. A clear evolution of 

the complexity of the field is present. The signed measure of the current can be then 

computed as usual: 



fii(r) = I J z (x,y) dxdy , 



and the scaling properties of the time averaged partition function are reported in figure [| figure [5| 

The power-law scaling (^) is clearly visible in a range extending from the large scales 

(near the integral scale of the flow £ ~ 0.2L, L = 2tt being the size of the simulation 

box) down to a correlation lenght r* of the order of the Taylor microscale A ~ 0.02L of 

the flow (see for example Frish, 1995). In this region, we fit the partition function to 

obtain the cancellation exponent k = 0.43 ± 0.06. A saturation of the partition function 

is observed at a scale r$ which is found to be of the order of the dissipative scale of the 

flow. In fact, for scales smaller than r$ the dissipation stops the structures formation 
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cascade, so that cancellations are stopped too. The fractal dimension of the current 
structures has been computed using the relation (H), which gives D ~ 1, indicating 
current structures similar to filaments. The presence of filaments can be clearly observed 
by a direct inspection of the current field contour plot, confirming the reliability of the 
model (see Sorriso-Valvo et ah, 2002). 

In order to compare more directely the numerical results with the solar data 
cancellation analisys, we introduce now a new surrogate for the current helicity. Since in 
two-dimensional geometry the current helicity is zero, the current being perpendicular 
to the magnetic field, we simply consider Hf? d \x,y) = J(x,y)\B(x,y)\, which is 



represented in Figure ^ for the same times snapshots as in Figure f|. The current helicity Figure ^| 
surrogate, as in the case of the solar data, appears smoother than the current itself, and 
keeps the same structures topology and evolution. The signed measure of such field is 
then computed as in the previous cases: 

Hi{r) = [ Hf d \x,y)dxdy. 

JQi{r) 

The scaling properties of the partition function x{ r ) can be now represented by the 
cancellation exponent, obtained by the usual fitting procedure after time averaging. 



In Figure [7| we present the scaling of x{ r )i together with the fitting power-law, for which Figure [7] 
we find an exponent k = 0.46 ± 0.03. This result is very close to the result obtained 
for the current, showing that for our numerical data as well, the current field is the one 
responsible for sign singularities, and its structures control the cancellations. 

We want now to study in more detail the time evolution of the cancellations effects 
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in the two-dimensional numerical simulations. To do this, we plot in Figure |8| the 
time evolution of the (kinetic) Reynolds number, together with the two cancellation 
exponents Kj and K H (2d) , for the snapshots already presented in Figures |] and As 
already pointed out, the first snapshots looks smoother than the others (see Figures f| 
and H and this fact can be interpreted as stronger presence of dissipative effects for such 
times; correspondingly, the values of both the cancellation exponents are smaller, which 
means, following our model, that the fractal dimension of the structures D is larger 
for these times. This observation confirm the results already obtained in Sorriso-Valvo 
et al, (2002) and discussed above. The Reynolds number also presents a time evolution, 
which is however shifted with respect to the evolution of the cancellation exponents. 
Unfortunately, due to the limited time interval of our simulations, it is impossible to say 
whether that shift is backward or foreward. Since the Reynolds number is related to the 
importance of dissipative effects against non linear effects in the turbulent cascade, it 
would be interesting to clarify this point as a further confirmation of our interpretation. 

In this paper we point out that the changes in the scaling behavior of cancellations, 
measured through the cancellation exponent k, are due to the topology changes of the 
structures present in the field, and are thus related to the importance of dissipative 
effects. The non-linear turbulent cascade, underlying the formation os such structures 
on all scales, can be considered as one important input mechanism for flares. The 
results obtained from the analysis of the numerical simulations strongly support the 
interpretation of the observational results for the photospheric magnetic field in the 
active regions. 
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To conclude, it is evident that the behavior we found can be used as a signature 
of the occurrence of big flares. High energy solar flares become of great interest 
because they can produce severe damages on Earth. Power blackouts, break up of 
communications and mainly damage of satellites or space flights, can be ascribed 
to energy released during big solar flares. It is then evident that the possibility of 
forecasting, even if partially, high energy flares has a wide practical interest to prevent 
the effects of flares on Earth and its environment. We build up a model which allows 
us to recognize without ambiguity changing behavior of the photospheric magnetic 
field of active regions. These changes, pointed out through the variation of a scaling 
index for current helicity, can be seen mainly before the eruption of big flares. The 
change of scaling index is due to the turbulent and intermittent energy cascade towards 
smaller scales, a mechanism which could be identified as the input of flaring activity, 
where energy is dissipated. The method could allows us to forecast, in real time, the 
appearence of the strongest flaring activity above active regions. 
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Figure 1. The current helicity H c measured for the active region NOAA 7590. The 
five different plots are taken at different times, and include a 50 hours interval. The 
presence of positive and negative strucutres on all scales is evidenced. The flat portion 
of field near the corners are due to removing in the projection effects. Over each plot the 
measured cancellation exponent k and the fractal dimension of structures D are reported 
(see further for explaination). 
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Figure 2. The scaling of the partition function for a flaring active regions (NOAA 7315), 
which started to flare on October 22, 1992. The power-law fit is indicated as a dotted 



line. 
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Figure 3. We present one flaring event observed in 1993 (see text for description) in the 
bottom part of the plot (line bars, in arbitrary units). The corresponding time variation 
of the fractal dimension D is reported (symbols). 
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Figure 4. The current field J obtained from high resolution two-dimensional numerical 
simulation of MHD equations. The different plots are ten snapshots in the statistically 
steady state, from about t = 168 up to t = 336 in non-linear times units, r NL . As for 
the solar data, the presence of positive and negative strucutres on all scales is clear. We 
report, over each plot, the measured values of k and D. 



20 



X(r) 10 ' 
10 2 





000 00 °^ 


/c=-0.43 




- 




- 






u 




























c 



10 3 10 2 10 1 10° 

r/L 

Figure 5. The scaling of the partition function for the current obtained from numerical 
data. This result is obtained by averaging the time evolution, in order to increase the 
statistics. The power-law fit is indicated as a straight line. The scales are normalized to 
the simulation box size L = 2-7T. 
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Figure 6. The surrogate of the current helicity Hf?^ obtained from high resolution 
two-dimensional numerical simulation of MHD equations. The different plots refer to the 
same snapshots presented for the current in Figure §. The signed structures are present 
and reproduce the current structures, but the field looks smoother than the current field 
itself. Again, we report, over each plot, the measured values of k and D. 
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Figure 7. The scaling of the partition function for the current helicity obtained from 
numerical data. As for the current, this result is obtained by averaging the time evolution, 
in order to increase the statistics. The power-law fit is indicated as a straight line. The 



scales are normalized to the simulation box size L = 2tt. 
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Figure 8. The cancellation exponents k for both the current (black circles) and the 
current helicity (stars) for different times. The open circles plot represents the kinetic 
Reynolds number of the flow (left y-axis). 



